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We consider a massless, minimally coupled scalar with a quartic self-interaction 
which is released in Bunch-Davies vacuum in the locally de Sitter background of an 
inflating universe. It was shown, in this system, that quantum effects can induce 
a temporary phase of super-acceleration, causing a violation of the Weak Energy 
Condition on cosmological scales. In this paper, we investigate the system's stabil- 
ity by studying the behavior of linearized perturbations in the quantum-corrected 
effective field equation at one- and two-loop order. We show that the amplitude of 
the quantum-corrected mode function is reduced in time, starting from its initial 
classical (Bunch-Davies) value. This implies that the linear perturbations do not 
grow, hence the model is stable. The decrease in the amplitude is in agreement with 
the system developing a positive (growing) mass-squared due to quantum processes. 
The induced mass, however, remains perturbatively small and does not go tachyonic. 
This ensures the stability. 
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I. 



INTRODUCTION 



Present cosmological observations [1] do not exclude the possibility of an evolving dark 
energy equation of state w = p/p whose current value is less than minus one [2], i.e., a 
phase of super-acceleration. Although the data are consistent with w — — 1, which can be 
explained by a simple cosmological constant, the possibility of w < — 1 has been an area of 
great interest in recent years [3]. 

Super-acceleration is difficult to explain with classical models on account of the problem 
with stability [4]. One can achieve models exhibiting w < —1, by postulating scalar fields, 
for example. Such models, however, decay irrespective of how this is achieved. The ob- 
served persistence of the universe, therefore, can only be consistent with a relatively brief 
self-limiting phase of super-acceleration. One way to get such a self-limiting phase, with- 
out violating classical stability, is via quantum effects [5-11]. The energy-time uncertainty 
principle requires virtual particles to emerge from the vacuum and then disappear back into 
it. The inflationary expansion of spacetime, however, causes the virtual particles persist 
longer than the flat spacetime [12]. In fact, any sufficiently long wavelength virtual particle- 
antiparticle pairs, which are massless on the Hubble scale, are pulled apart by the Hubble 
flow before they find time to annihilate each other. Hence, they become real and may persist 
forever; recalling the analogy with the Hawking radiation. The rate at which the virtual 
particles emerge from the vacuum, on the other hand, is suppressed by the inverse of the 
scale factor for conformally invariant particles. Thus, quantum effects are enhanced during 
inflation for particles that are effectively massless (with respect to the Hubble parameter H) 
and classically conformally noninvariant. Gravitons and massless minimally coupled (MMC) 
scalars are unique in possessing zero mass without having classical conformal invariance. One 
is lead, naturally, to a self-limiting quantum effect in a classically stable theory, such as the 
MMC scalar with a quartic self-interaction in the locally de Sitter background of an inflating 
universe. The Lagrangian density that describes this system is 



The dynamical variable in the model is the scalar field <p(x). The metric g^ v is a nondy- 
namical background which is taken to be a D-dimensional locally de Sitter geometry. The 



C = --^-gg^d^d uV 




gip 4 + counterterms . 



(1) 
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invariant element can be expressed conveniently either in comoving or conformal coordinates 

ds 2 = -dt 2 + e 2Ht dx- dx = a 2 (r]) -drf + dx ■ dx , (2) 

respectively. The conformal factor and the transformation which relate the two coordinate 
systems are 

a( V ) =-±. = . (3) 

The Hubble constant H is related to the cosmological constant A = (D — 1)H 2 . It is the 
cosmological constant that drives inflation in the model. The scalar is a spectator to A-driven 
(de Sitter) inflation. We adopt the following notations: x 11 = (x°,x), x° = rj, = (<9 , V). 

We release the state in Bunch-Davies vacuum at t — 0, corresponding to conformal time 
i] = rji = —H~ l . Hence, the scale factor is normalized to a = 1 when the state is released so 
that a > 1 throughout the evolution. Note that the infinite future corresponds to rj — > 0~, 
so the possible variation of causally related conformal coordinates in either space or time 
is at most Ax = Arj = H' 1 . Applying the Schwinger-Keldish formalism [13, 14] and using 
dimensional regularization, the fully renormalized vacuum expectation value (VEV) of the 
stress-energy tensor, is calculated [8, 9] in this system. The energy density 

p = (Q\T 00 \Q) / a 2 (rj) and pressure p5y = (Q\Tij\Q) / a 2 (rj) are obtained as 

a- - ^ + ^{5 h >> + l«"-5|:<^«^ +0 <*>- < 4 » 



A \H 4 l... 1 . . l^n 2 -4 (n+1) ) ^ /a2 , 

p = —{ - In 2 (a) + - In (a) + - > -=a ( n+ 1] \ + 0{\ 2 ) 

y 8nG 2 6 tt 4 1 2 v ; 3 V ; 6 ^ (n + 1 2 j K ' 

v n=l v ' ' 



(5) 



Notice that p rcn and p Icn obey [8, 9] the covariant conservation law T lJLV {x) ' ,u = 0, i.e., p rcn = 
— 3H(p rcn +p Icn ), where the dot denotes derivative with respect to the comoving time t. Their 
sum, however, violates the Weak Energy Condition (WEC) p+p>0 on cosmological scales 



XRA f 1 , / \ 2 _ 3 l^n + 2 , n+1) \ 

^ n=l J 



)• (6) 



Although the value for w + 1 is unobservably small in this model, the calculation shows 
that quantum effects can induce a self-limiting phase of super-acceleration in which a clas- 
sically stable theory violates the WEC on cosmological scales in the average of p+p, not 
just in fluctuations about an average that obeys the condition p + p = 0. This is be- 
cause inflationary particle production causes the scalar to undergo a random walk such 
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that its average distance from the minimum of the potential ^p 4 increases. In our model, 
(n\ip 2 (x)\n) = (UV divergence) + H 2 ln(a)/47r 2 + 0(A) [15]; recall that \n(a)=Ht. (See the 
calculations in Sec. IV for the 0(A) and 0(A 2 ) corrections.) Hence, after the ultraviolet 
divergence is removed, the VEV of p 2 gets pushed up its potential by inflationary particle 
production. This increases the vacuum energy which leads to the violation of the WEC by 
virtue of the covariant conservation p = —3H(p + p); since p > due to inflationary particle 
production, p + p has to be less than zero. 

The process, however, must be self-limiting because (i) as the scalar rises up its potential, 
the classical restoring force — A</? 3 /6 pushes it back down, and (ii) the curvature \p 2 /2 
associated with being away from the minimum of the potential acts like a positive "mass- 
squared" to reduce the inflationary particle production responsible for pushing the scalar 
away from the configuration tp — where the potential is minimum. (In quantum field 
theory (QFT) the mass-squared is calculated via self-energy diagrams as is rigourously done 
in Ref. [10]. The VEV of the curvature of the potential provides a heuristic picture to 
understand the effect.) Since the classical restoring force (i) gets bigger as the field rolls 
up its potential and the mass generation (ii) cuts off particle production, the field cannot 
continue rolling up its potential. It must eventually come to a halt. Indeed, Starobinsky and 
Yokoyama showed [7] that (Q\p 2 (x)\Q) asymptotes to the constant 3H 2 T (3/4) /itT(l/A)\/\ 
in this model, which proves that the field strength does not grow forever. The curvature of 
the potential, that acts like mass-squared, should asymptote to A/2 times this expectation 
value. They also estimated the time scale for the process as T xs 18.7/ Hy/X. Thus, by 
choosing A <C 1, it is possible to have long duration for the effect. We assume A 1 in this 
paper. 

We study the stability of the system in this paper. To decide whether the system is stable 
[15] or not, one needs to check (i) if the VEV {£l\(p 2 (x)\il) continues to grow without a bound, 
and (ii) if the small, position-dependent perturbations grow. If neither happens, the system 
is stable; otherwise, it is unstable. The above arguments show that (£l\p 2 (x) \Q) cannot 
continue to grow forever in the interacting theory (it asymptotes to a constant). Checking 
criterion (ii) is the main object of this paper. To do that, one solves the quantum-corrected 
effective field equation at linearized order 




(7) 
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and obtains the quantum-corrected mode function. Although the scalar is classically mass- 
less in our model, quantum processes generate a nonzero self-mass-squared M 2 (x;x'). Po- 
tential instabilities would come from the field developing a negative mass-squared. In that 
case, the amplitude of the mode function would be an increasing function of time, indicating 
growth of perturbations and, hence, the instability. The fully renormalized scalar self-mass- 
squared M 2 {x\ x') is calculated rigorously in Ref. [10] at one- and two-loop order, using the 
Schwinger-Keldish formalism. M 2 (x;x') is indeed positive at one- loop. However, one must 
go to two-loop order to see corrections of the derivative terms. To interpret the two-loop 
result, and hence to check the stability of the system, one needs to investigate how the 
self-mass-squared M 2 (x;x') modifies the effective field equations and its solution, i.e., the 
quantum-corrected mode function. If the amplitude of the solution is a decreasing function 
of time, one can conclude that perturbations do not grow; therefore, the model is stable. 

The outline is as follows. In Sec. II, we define the effective mode equation, summarize 
the Schwinger-Keldish formalism, and discuss our limitations in solving the effective mode 
equation. In Sec. Ill, we solve the effective mode equation in late time limit and obtain the 
mode function in leading logarithm approximation. Late time, for us, means ln(a) ^> 1. In 
Sec. IV, we alternatively compute the same mode function using Starobinsky's stochastic 
inflation technique, and compare it with the result obtained in Sec. III. Our conclusions are 
summarized in Sec. V. 

II. EFFECTIVE MODE EQUATION FOR THE MMC SCALAR 

In this section, we describe the operator formalism and effective field equation correspon- 
dence. Then, we review the Schwinger-Keldish formalism that one must use to calculate 
expectation values. We use the one- and two-loop results [10] for the scalar self-mass- 
squared M 2 (x;x'), obtained by applying the Schwinger-Keldish formalism in our model, to 
write down the effective (quantum-corrected) mode equation and discuss how we "solve" it. 

A. Relation to Fundamental Operators 

The relation between the fundamental Heisenberg operator of scalar field (f(x), and the 
C-number plane wave mode solution Q(x; k) of the linearized effective field equation can be 
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given [16, 17] as 



k) 



(8) 



Here and are the states, and cr(&;) is the free creation operator. In flat space 
scattering problems, and correspond to the states whose wave functionals are free 
vacuum in the asymptotic past and future, respectively. The universe, however, begins at 
a finite time, and evolves to some unknown state in the asymptotic future. Therefore, in 
cosmology, we release the universe from a prepared state at a given time and then let it 
evolve. We seek to know expectation values in the presence of this state. This corresponds 
to the choice = For computational convenience, we assume that both of the 

states are free vacuum at r] — r)i. Had the choices of flat space scattering theory been 
used, acausal effective field equations would have been obtained. The matrix elements of 
Hermitian operators would also be complex in that case. 

To define the free creation and annihilation operators, recall that the full Lagrangian 
density C of the MMC scalar is 



(1 + 5Z) 



(A + 5X) 



5m? 



~9<P 



(9) 



The field strength (5Z), coupling constant (SX), and mass (5m 2 ) counterterms are needed 
to remove divergences at one- and two-loop order in the scalar self-mass squared. It turns 
out that 5Z and 5\ are of order A 2 , whereas 5m? has contributions of order A and A 2 [10]. 
Let us now integrate the invariant field equation of the MMC scalar 



-9 



+5Z 



(X + 5X) 3 . 2 
6 



= 



(10) 



The result is 



(p(x) = M*)+ [ drf ! d D ' 1 x'G(x]x')I[ip(x') 



(11) 



where ^q(x) is the free field. We define the interaction term as 



+5Z 



(X + 5X) o P , 
-(f 6 + 5m 2 (p 



(12) 



The Green's function G(x; x') is any solution of the equation 

d,(y^g^d v G(x ]X ')) = 5 D (x-x') . 



(13) 
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Although the Green's functions would obey Feynman boundary conditions for flat space 
scattering problems, it is more natural to use retarded boundary conditions in cosmology. 
The fundamental field operator (p(x), on the other hand, is unique. It does not depend on 
the choices of the boundary conditions for the Green's functions or on rji. What changes 
with those choices is the free scalar ipo(x). Because <^o(x) obeys the linearized equations of 
motion and agrees with the full fields at rj — rji, it can be expanded in terms of free creation 

— » — * 

and annihilation operators a(k) and a*(k) as 

MX) = J (CpM^' fc ) e ** a $ + *)e-*Mft} , (14) 
where the Bunch-Davies mode function [18] 

u(r),k) = (l +ikv)e~ ikv ■ (15) 

Although the creation and annihilation operators change as different Green's functions are 
used in Eq. (11), their nonzero commutation relation remains fixed 

[a(k),a*(k')] = (2n) D - 1 5 D - 1 (k-k') . (16) 

By iterating Eq. (11), one can expand the full field (p(x) in terms of the free field ipo(x) as 

<p(x) = <p (x) + [°H [ d D - l x'G^{x-x')I[ip{x')\ (17) 



= <p (x) + f drj' j d D ~ l x'G^{x- x')I[<p (x')} + ... . (18) 

Hence, choosing |^) = as free vacuum at rji, one can see [16, 17] that the quantum- 
corrected plane wave mode solution (8) yields 

$(x; k) = (tt (p(x), a\k) = u(r), fc)e* 2 + 0(A) . (19) 

The 0(A) and 0(A 2 ) corrections in Eq. (19) are obtained in Sec. Ill, by "solving" the 
quantum-corrected effective field equation at one- and two-loop order. In Sec. IV, we obtain 
the same corrections by calculating commutator (8) stochastically. The results yielded by 
the two approaches are in perfect agreement. 



B. Schwinger-Keldish Formalism 

Because of the fact that "in" (t — > — oo) vacuum is not equal to the "out" (t — > oo) 
vacuum in de Sitter background, we need to calculate expectation values, rather than in- 
out matrix elements. This is done by applying the Schwinger-Keldish formalism [13, 14]. 
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The endpoints of propagators acquire a ± polarity, in this formalism. Hence, every propa- 
gator iA(x;x') of the in-out formalism generalizes to four Schwinger-Keldysh propagators: 

iA ++ (x; x f ), iA + _(x; x'), iA_ + (x; x') and iA (x; x'). Each propagator can be obtained from 

the Feynman propagator by replacing the de Sitter conformal coordinate interval 

2 / \ 2 



Ax (x,x') = Ax ++ (x;x' 



with the appropriate coordinate interval, 



7/- 



(20) 



Ax 2 + _(x; x 1 ) 



x—x 



(rj-rj'+iSy = (Ax 2 _ + (x;x'))* , 



Ax 2 __(x;x') = (Ax 2 ++ (x;x'))* . 



(21) 
(22) 



Vertices are either all + or all — . A + vertex is the usual one of the in-out formalism, 
whereas the — vertex is its conjugate. 

Because each external line can be either + or — in the Schwinger-Keldysh formalism, 
each iV-point 1PI function of the in-out formalism corresponds to 2 N Schwinger-Keldysh 
iV-point 1PI functions. The Schwinger-Keldysh effective action is the generating functional 
of these 1PI functions, so it depends upon two background fields (p+(x) and <p-(x). For 
example, there are four Schwinger-Keldysh 2-point 1PI functions Ml ± (x;x'). The +- 1- one 
is the same as the in-out self-mass-squared and the others are related as the propagators 

-iM 2 _(x;x') = {-iM 2 ++ (x;x')y , -iM 2 + (x;x') = (-iMl_(x; x')^* . (23) 

The various self-mass-squared terms enter [14] the effective action as follows: 



F[<P+,<P-] = >%+] - S[p. 
_I J d D x J d D x' 



+ (24) 



^ + (x)M 2 ++ (x; x')(fi + (x') + ^ + (x)Ml_(x; x')^_(x') 
+ip^(x)M 2 + (x; x')ip+(x') + ip-(x)M 2 _(x;x')ip„(x') 

where S[ip] is the classical scalar action. The effective field equations of the Schwinger- 
Keldysh formalism are obtained by varying with respect to either polarity, and then setting 
the two polarities equal [14]. Up to order 0((f 2 ), we have 

5T[<p ± ] < ~ '° 

Note that we have taken the regularization parameter D to its unregulated value of D = 4, 
in view of the fact that the self-mass-squared is assumed to be fully renormalized. It is this 



dviy^gg^dM*))- f drj' f d 3 x'[M 2 ++ (x;x') + M 2 + _(x;x')}ip(x') . (25) 
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linearized effective field equation which k) (Eq. (8)) obeys [16, 17] 

dr/ y rfV|M^ + ( a ;; a ; , ) + M^(x;x , )}$(x , ;fc) =0. (26) 

Thus, the two renormalized 1PI 2-point functions we need are M 2 + (x;x') and M 2 _(x;x'). 
At one-loop order, we have [10] 



M 2 (x-x') 



XH 2 

8tt 2 



a 4 \n(a)5 4 (x - x) + 0(A 2 



(27) 



The H — case vanishes at this order because there is no mixed interaction. Fully renormalized 
two-loop results for the ++ and H — cases are [10] 



M, 



i\ 2 \ aa' 



2++ 



d 4 



In (/j, 2 Ax 2 ++ ) 



Ax 2 



-H 2 (aa') 2 d 2 



-H 4 (aa') 
X 2 



2% 6 \ 24 

ln 2 (^ 2 A, 2 ++ ) 
-3 __ ^ -(^(^A, 2 

— a 2 1 - ln(a)<9 2 + (2 ln(a) + i)#a<9 j^x-a;') 



(Heh In (/i 2 Ax 2 + ) In 2 (/i 2 Ax 2 + ) 



/tie*\ 



Ax 2 ++ 



4Ax 2 ++ 



+ 



2 9 3 

A 2 # 2 f 4, 3 . . 23. 2/ , 
2V -g ln(a) -18 ln(a) + 
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+ 3 In 



(|)-H in(a) } 

7 -("+2) 



X 2 H 2 f a" 3 



ln(a) |>a 4 5 4 (x-x') 

00 -(n+3) 



2 7 7T 4 



f r> 4- ^ ^ n- n + 2 n- n + 3 1 

1 81 ^ (n + 1) 3 ^ n + 2 3 ^n n + 3 3 J v ; v ; 

v n=l v ' n=l v ' n=l v ' ' 



M, 



2+- 



iX 2 { aa' 
'2%e\ 24 



9 4 



In [n 2 Ax 2 + _) 



Ax 2 



-H 2 (aa') 2 d 2 



^ He'iynj^Ax 2 ^) | In 2 (/i 2 A:r 2 _) 



V 2/j 



Ax 2 



4Ax 2 



-H 4 (aa'f 



In 2 f 4H 2 Ax 2 



\ 4 



Ax 2 



+ ^_( aa ')4 ln 3^ iJ 2 A;r 2 

6 V 4 



(29) 



The ++ and H — terms in (26) exactly cancel for rj' > r\ and also, in the limit 5^0, for 
x'^ outside the light-cone of x M . This is how the Schwinger-Keldysh formalism gives causal 
effective field equations. In the next section, we discuss what we mean by "solving" the 
quantum-corrected effective mode equation (26). The one- and two-loop corrected mode 
solution is obtained in late time limit, i.e., for ln(a) 1, in Sec. III. 



C. Solving the quantum-corrected effective mode equation 



Here we discuss the limitations that one has in solving the effective mode equation (26). 
The full scalar self-mass-squared can be expressed, as a series, in powers of the loop counting 
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parameter A 

oo 

M 2 ++ {x; x 1 ) + M 2 + _{x; x 1 ) = ^M]{x- x') . (30) 

The first limitation is that we have only the £=1 and 1 = 1 terms 

M 2 ++ (x;x') = XMj(x;x') , (31) 
M 2 ++ (x;x') + M 2 + _(x-x') = X 2 Ml(x;x') , (32) 

which are given by Eqs. (27), (28) and (29), respectively. So we can only solve the effective 
mode equation perturbatively. We first substitute a series solution of the form 

oo 

k) = u(r], k)e 1 ^ + \%(r}, k)e ils (33) 
e=i 

into Eq. (26) and then solve the equation order by order in powers of A and A 2 . The zeroth 
order (£ = 0) solution of $£ is the well-known Bunch-Davies mode function u(rj, k) (Eq. (15)) 
times the exponential e lk ' x . 

The second limitation is due to the lower bound "r^" on the temporal integration in 
Eq. (26). We release the universe in free vacuum at time rj = rji. Little is known about 
the wave functionals of interacting QFTs in curved space, but free vacuum can hardly be 
realistic. In fact, all of the finite energy states of interacting flat space QFTs have important 
corrections. Similar corrections are expected in curved space, too. Although it is possible 
to correct the free state functionals perturbatively as in nonrelativistic quantum mechanics, 
the usual procedure in flat space QFT is to release the system in free vacuum at asymptotic 
past, and let the infinite time evolution resolve the difference between free vacuum and 
true vacuum into shifts of the mass, field strength and background field [19]. In cosmology, 
however, one cannot typically apply this procedure, for the reasons noted in Sec. II A. 
One can still correct the state wave functionals perturbatively, though. Corrections to the 
initial state would appear as new interaction vertices on the initial value surface. They are 
expected to have a large effect on the expectation values of operators near the initial value 
which would decay in the late time limit. For example, it is the exponentially falling portions 
of the renormalized stress-energy tensor (4) and (5) 

A,,^^{^4E^«-*" +1, } + ^). (34) 

v n=l v ' ' 

^ = ^{-l±^- ( " +1, }+0^ (35) 

n=l ' 
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that it is conjectured [9] can be absorbed into an order A correction of the initial (a = 1) 
free Bunch-Davies vacuum state. The fact that they fall off as one evolves away from the 
initial value surface suggests that they can be absorbed into a kind of local interaction there, 
leaving only the infrared logarithms 

^ = 8^ + ^{l' n2 < a) } +0(A2)l <36) 

Notice that they are separately conserved, i.e., p conj = —3H(p conj + p conj ) and p famng = 
— 3if (p falling + Pfamng). This is exactly what would be the case if they could be cancelled 
by a new interaction vertex. Note that Eqs. (34) and (35) diverge on the initial value sur- 
face at a = 1 which indicates that free vacuum is very far away from any physically accessible 
state. Thus, although Eq. (26) determines the quantum corrections to the mode function 
(19) for free vacuum, that mode function has little physical relevance, because free vacuum 
is inaccessible. To find physically relevant mode functions, which are also valid for initial 
times, the corrections to the state wave functional must be included. Unfortunately, we have 
neither order A nor order A 2 corrections to the state wave functional. It therefore makes no 
sense to solve Eq. (26) for all times. The effects of the state corrections, however, must fall 
off at late times (ln(a) 3> 1) as in Eqs. (34) and (35). Because of time evolution, initially 
free vacuum and true vacuum become indistinguishable, as in flat space QFT [19]. Hence, 
we may obtain valid information from Eq. (26) by solving it in late time limit. That is the 
subject of the next section. 



III. EFFECTIVE MODE FUNCTION FOR THE MMC SCALAR 

The linearized effective field equation that the MMC scalar mode solution $(x, k) obeys 
is given in Eq. (26). Using Eq. (30), (31) and (33), one obtains the integro-differential 

— * 

equation for the one-loop correction $1(77, k) to the classical mode function 7/(77, k) 
a 2 [d 2 + 2Had + k 2 }^ 1 (r ] , fc) = - I drf I <Px'M\{x- x')u(r)' , k)e~ in < 3 -^ 

Jr)i J 

f° r h 2 - h 2 

= - drf d 3 x'— a A \n{a)5\x-x')u{7 1 ' ) k)e- ik - { ^' ) = -7/(77, k) a 4 ln(a) . (38) 

As is discussed in Sec. II C, the only sensible and physically interesting regime in which we 
can solve the effective mode equation is the late time limit ln(a) 3> 1. The zeroth order 
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mode function u(r], k) can be replaced by its limit u(0, k) = H/ \Jlk z in this regime. Solving 
Eq. (38) in late time limit, we find 



(39) 



in leading logarithm orders. 

The order A 2 correction $ 2 (^, on the other hand, has contributions due to both one- 
and two-loop self-mass-squared terms. It obeys 



a 2 [d 2 + 2Had + k 2 ]^ 2 (r],k) 

= - f drf f d 3 x' {M 2 1 (x;x')^ 1 (r ] ',k)+Ml(x;x')u(r] , ,k)} 

Jr)i J 



-ik-(x-x') 



(40) 



The first integral is evaluated, in leading logarithm order, by inserting Eqs. (27) and (39) 
into Eq. (40). We find 



f drf I dPx'Mlix-x')®^'^) 



H 2 



2 7 3tt 4 



u(0, k)a 



ln 3 (a) - -ln 2 (a) 



(41) 



The second integral is evaluated, in late time limit, in Appendix A. Expanding in terms of 
powers of infrared logarithms, we find 

-0 r r0 



- I drf I d 3 x'M 2 2 (x;x')u(r]',k)e- tH£ - £,) — > -u(0,k) \ drf \ d 3 x'M 2 2 (x;x') 

23 



H 2 

~ — — = — tu(0, k)a 4 



ln 3 (a) + - ln>)+ - ln (| 



2 4 3 2 7T 4 

Using Eqs. (41) and (42) in Eq. (40) yields 
a 2 [d 2 + 2Had + k 2 ]<5> 2 (r], k) 



27, /2u\ 189 ?H ,. 

1 ' ' ' ^2- + y) ln(a) 



(42) 



H 2 



-w(0, k)a 4 



2 7 3 2 7T 4 

In leading orders, the solution for $ 2 (r], k) is 
1 , fl 1 



11 ln 3 (a) + f ln 2 (a) + (271n(|) - ^ + 4vr 2 ) ln(a) 



(43) 



$2(77, k). 



2 8 3 3 vr 



-u(0,k) <Jy ln 4 (a)-ln 3 (a) + 



27 In 



2/i\ 185 



(?) 



f47T / 



ln 2 (a)--ln(a) 



.(44) 



Thus, keeping the leading logarithm terms at each order of perturbation (i.e., in <!>! and $2)) 
we find that the quantum-corrected mode solution (33) asymptotes to 



<3>(x, k) ~ u(0, k)e l 



1 - 



2 4 3vr : 



: Aln 2 (a) + 



11 



2 9 3 3 7T 4 



A 2 ln 4 (a) \ + 0(A 3 ) . (45) 
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One can immediately see from Eq. (45) that perturbation theory breaks down when In (a(i)) 
is of order 1 / y/\. This, however, does not invalidate the reliability of our late time (ln(a) 3> 1) 
solution, because by choosing A <C 1, as we assume in this paper, one can have a long period 
of time during which 1 <C ln(a) <C 1/VA. 

Equation (45) also shows that, at t — 0, the mode solution k) is equal to the well- 
known (Bunch-Davies) classical result u(0, k)e lk ' x . As time goes on, it decreases proportional 
to the factor 1 — \H 2 t 2 /A8n 2 + 0(A 2 ) (the stochastic calculation of Sec. IV yields the same 
result). Thus, the amplitude of quantum-corrected mode function (hence, of the field) is 
reduced (consistent with the model developing a positive mass-squared, as one- and two- 
loop self-mass-squared terms and the VEV of the curvature of the potential imply; see the 
discussion in the next section). This means that linear perturbations do not grow in this 
system; therefore, it is stable. 



IV. STOCHASTIC ANALYSIS 

Starobinsky developed a stochastic inflation technique [6, 7] which gives the leading 
infrared logarithms at each order in perturbation theory. Recently, his technique was proven 
to all orders and extended to various models [20]. In this section, we introduce the stochastic 
technique briefly, and use it to calculate the quantum-corrected plane wave mode solution (8) 
and the VEV of the curvature (associated with the field being away from the minimum) of 
the potential which acts like a mass-squared in the classical action. 

The equation of motion for the scalar field with quartic self-interaction in D = 3 + 1 
dimensional de Sitter background is 

V 2 A 

Cp{t,x) + 3H(p(t,x) - -r<p(t,x) + -y? 3 (t,x) = . (46) 

cr 6 

The solution of Eq. (46) can be obtained by iterating 

<p(t, x) = <p (t, $)-^f <*fa 3 (0 J d 3 x'G iet (t, x- 1', x') V 3 (t', x') , (47) 
where the retarded Green's function [20] is 

^ 5 ^ t - [ j(g!!^ff^ +e(g || g -^|| +B -. (t) - a -. (0) ] . (48) 
As in Eq. (14), the free field (po(t, x) can be expanded in terms of mode function u(t, k) and 

— » — * 

annihilation and creation operators a(k) and a'(k), satisfying the canonical commutation 
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relation (16) in D = 3 + 1 dimensions. Starobinsky's stochastic technique cuts out the 
ultraviolet modes k > Ha of the field and applies the following rules to the equation of 
motion: (i) retain only the term with the smallest number of derivatives of the field, (ii) 
replace the field variable by a stochastic variable, and (iii) subtract the stochastic source 
term / for each time derivative of the field. Applying these rules to Eq. (46) yields 

(p + 3H<p - ^ + \<P* = — SHip + ^ 3 = — 3H (0 - / ) + = . (49) 

(The scalar field cp became a stochastic field 4>.) The source term is the time derivative 
of the infrared truncated free field (14) 

Mt,x) = J-0jS - k ) ;fp { e &3 a(k) + e~ &3 a\k) } • (50) 

Here the leading infrared limit of the Bunch-Davies mode function (15) u(t,k) ~ H/\ / 2k 3 
is used. Hence, 

U(t,x) = Mt,x) = J—5{k-Ha{t))^^ k -*a{k)+e-^a\k)y (51) 
In Eq. (49), we obtained a Langevin-like equation which can be recast as 

0(t,f)=^(t,f)- I ^:0 3 (t,f). (52) 

In this section, we use this equation to stochastically calculate (i) the VEV of the curvature 
associated with being away from the minimum of the potential, i.e., \(Q\(fi 2 (x)\Q) /2, which 
acts like field-dependent mass-squared, and (ii) the quantum-corrected mode function (8), 

— * 

i.e., (Q\[(f)(x), a*(k)]\Q). As a check, we calculate (i) also using perturbative QFT at one- 
and two-loop order and show that the two realizations agree perfectly in leading logarithm 
order. The quantum-corrected mode function was already obtained in Sec. Ill by applying 
QFT. Comparing Eq. (45) of Sec. Ill with the stochastic result for (ii) will show that the 
agreement is again perfect in leading logarithm order. 

In stochastic calculations (i) and (ii), we express (f>(x) in terms of the infrared truncated 
free field 0o perturbatively, by first integrating Eq. (52) and then iterating the result up to 
the desired power of A: 

0(t, x) = Mt, x ) --^ jf dt'^{t\ x) 

= Mt, x) - J*dt><f>l(t', x) + Jdt'^it', x) J*dt"4(t", x ) + 0(A 3 ) . (53) 
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In the noninteracting (free) theory, (j>(t,x) = (f) (t,x). Hence, using Eq. (50), the VEV of 
the scalar field-strength-squared is obtained trivially in this (A = 0) limit 



H 2 



<n|^(x)|n> = i ^in(a). 



(54) 



This stochastic result is the same as the results of Refs. [15] applying QFT. 

Now we start calculating (i) the VEV \(Vt\4> 2 (x)\Vl) /2 in the interacting (A ^ 0) theory. 
Using Eq. (53) we find 



~ 2 



(n 1 0o (x) | n> - Ajf W (n | ^ (t, x) 4 (t ,x)\n) 



+0(A 3 ) 



(n|^g(x)|n>- 



A 



3H 



dt'(n\<f) (t, x)(f) (t', x)\n) (n\(f)l(t', x)\n) 



+0(A 3 ) .(55) 



Inserting Eq. (50) into Eq. (55) yields the stochastic result 

H 2 



^|0 2 (:r)|fi> = ^ 



4tt 2 



ln(a) 



H 2 

2'V 



■A ln(a) 



3£T/ " 16 7T 4 

1 



ln 2 (a / )+0(A i ) 



2 2 3% 2 



Aln 2 (a) 



-0(A 3 



(56) 



Next, we want to calculate the same VEV A(O|0 2 |O) /2 using QFT. Figure 1 depicts the one- 
loop contribution. Hence, at one-loop order, the VEV is given in terms of the coincident 




8m 2 

® 



FIG. 1: Generic one- loop diagram with mass counterterm 



limit of the scalar propagator and the mass counterterm 5m 2 : 



^{n\(f) 2 \n} = ^iA(x; x) + 5m 2 + 0(A 2 ) . 



The scalar propagator in D-dimensional locally de Sitter background is [8, 9] 



' ] (4vr)#l ^n-f + ir( n + l) W 



iA(x: 



1 T(n+ f ) (yy-%+^ F(D - 1 



(57) 



D . 



^— TTCOt (7T — ) 



1 F(n + D- 1) r(D-l) , M l 

ir(!) +J F(ir ln ( aa) V (58) 



n=l v 2 
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Here the modified de Sitter length function y(x; x') is given in terms of the de Sitter conformal 
coordinate interval Ax 2 (Eq. (20)) 



y(x; x') = H 2 aa'Ax 2 = H 2 aa' [\\x - f || 2 - ( j 77 - rf\ - i5) 2 ] 



(59) 



To facilitate dimensional regularization, we express the dimension of spacetime in terms of 
its deviation from four: D = 4 — e. Therefore, the coincident limit of the scalar propagator 



iA(x;x) = lim iA(x;x') 



H 2 ~ e r(3-e) f rtl , , /vre 



(60) 



Because of the finite, time-dependent term in Eq. (60), we cannot make the one- loop 
VEV (57) vanish for all time. Our renormalization condition is that it should be zero 
at t — 0, which implies 

\H 2 ~' r(3-e) 



5m 2 = 

Therefore, Eq. (57) yields 



2 5 - e 7T 2 ^ 



r(2-§) 



-^ 7 rcot(^)+0(A 2 ). 



(61) 



(62) 



The two-loop diagram that contributes to the VEV is known as the snowman diagram 
depicted on the left of Fig. 2. The right-hand side diagram depicts the one-loop mass 




FIG. 2: Generic snowman diagram with mass counterterm. 

counterterm which naturally combines with it (5m 2 denotes mass counterterm vertex). In 
Schwinger-Keldish formalism (Sec. II B) the internal vertices are summed over both + and 
— polarities. A simple application of Feynman rules gives 



-iX) 



iA{x';x')-iSm 2 }. (63) 



^(fi|0 2 |fi) 2 _ loop = ^xV^A ++ (x;xO] 2 -[^A + 4x;xO] 2 
Both ++ and H — propagators are the same function (58) of the appropriate version of the 



modified de Sitter length function y(x;x'). By definition (20), y ++ (x;x') = H 2 aa'Ax 



++ 
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y(x;x'), given in Eq. (59). On the other hand, y + -(x;x') = H 2 aa'Ax 2 + _, where the coordi- 
nate interval Arr^_ is given in Eq. (21). The coincident propagator and the mass counterterm 
are calculated in Eqs. (60) and (61), respectively. Because both diagrams in Fig. 2 have the 
same lower loop, they posses the common factor given in the first curly bracket of Eq. (63). 
The first term in the second curly bracket comes from the left-hand side diagram, whereas 
the second term comes from the right-hand side diagram. 

The integral in Eq. (63) is calculated explicitly in Ref. [10]. The result can be read off 
directly from its Eq. (61). After renormalizing the overlapping divergence — ^r(jf^) el -^- 
of the snowman diagram by the two-loop mass counterterm, one obtains 



1 2 1 0) 2 _ loop = ^ (- 1 In 3 (a) + 1| In 2 (a) + 



2 7 tt 4 9 v ' 18 



VW 3 



, , \ 238 13 2 

ln(a) 1 it 2 

K J 81 54 



A „-3 00 „ i K 00 n -(n+2) 00 n -(n+3) ^ 

n=l x ' n=l v ' n=l ' 

where the Euler's constant 7 ~ 0.577. (The constants and exponentially decaying terms 
may be subsumed into the definition of the vacuum state.) Combining one- and two-loop 
results, Eqs. (62) and (64), we find 



= J^Aln(a) [l--^Aln 2 (a)j+0(A 3 ) , (65) 

in leading logarithm order. This QFT result is exactly the same as stochastic result Eq. (56). 

Until the breakdown of the perturbation theory — that occurs around A In 2 (a) ~ 1, as 
was pointed out earlier — expectation value (56) remains positive. This means that, the 
curvature associated with the scalar being away from the minimum of the potential assumes 
a growing positive expectation value which acts like a positive "mass-squared" during the 
process. This agrees with the decreasing mode function obtained in Eq. (45) by solving the 
one- and two-loop corrected effective field equation, in the context of QFT. 

Next, in calculation (ii), we recompute the very same mode solution $(x; k) obtained in 
Eq. (45) using QFT. This time, however, we calculate k) = (Q\[<j>(x), aJ(k)]\Q) (Eq. (8)) 
by applying the stochastic technique to show that the two realizations yield exactly the same 
result, in leading logarithm order. The commutator 

ft \2 



0(x), a\k)] =«(0, fc) e ** jl - A J <%(!?, + dt'<f>l(t', x) J dt"<t> 2 (t", x) 

jdt'4> {t\ x) J dt"(f)l(t", f)| + 0(A 3 ) , (66) 



A 2 

+ 2 3 3 H 2 
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which implies 



A 



dt'(n\(f)l(t',x)\n) 



A 2 



2 2 32^2 
t r t' 



dt' / dt" (n\<i) 2 {t',x)\tt)(n\<f) 2 (t" , x ) 1 n) 

Jo 



+2 / dt' df[{n\<j>o(t',x)<l>o(t n ,x)\n) 
Jo Jo 



+ 



A 2 



2 3 2 H 2 

Using Eq. (50) in Eq. (67), yields 



dt' dt''(Q\(Po(t',x)Mt''ix)\ty(n\(pl(t'',x)\n) \ +0(A 3 ) . (67) 
Jo 



£;)=u(0,£;)e**<{l 



A 

OH 



'dt'^Ha') 



A 2 

+ 2 2 3 2 H 2 

u(0, ky t3 



1 



,-Aln 2 (a) + 



11 



A 2 ln 4 (a) |> +0(A 3 ) 



(68) 



2 4 3tt 2 ' 2 9 3 3 tt 4 ' 

This result is in perfect agreement with Eq. (45), which is obtained by lengthy and highly 
nontrivial calculation using quantum field theory. 



V. CONCLUSIONS 

Massless, minimally coupled ^y? 4 on a locally de Sitter background can induce enhanced 
quantum effects causing super-accelerated phase of cosmic expansion, a possibility not ex- 
cluded by present observations. In this paper, we have studied the stability of this system 
for A <C 1. In Sec. II, we have obtained the quantum-corrected effective field equations at 
linearized order, using the fully renormalized Schwinger-Keldish self-mass-squared terms at 
one- and two-loop orders. In Sec. Ill, we have solved the effective field equations in the 
late time limit, i.e., for ln(a) ^> 1, and obtained the scalar mode function in leading powers 
of infrared logarithms at each order of perturbation. In Sec. IV, we have used Starobin- 
sky's stochastic inflation technique to compute the mode function in the leading logarithm 
approximation and compared it with the quantum field theory result of Sec. III. 

There are three main conclusions that we draw: (i) perturbation theory breaks down for 
ln(a(t)) ~ 1/VX This, however, does not invalidate the reliability of our late time solutions 
since one can have long period of time during which 1 <C ln(a) <C 1/a/A for A <C 1. (ii) The 
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quantum-corrected mode function decreases in time — consistent with the field developing 
a positive (nontachyonic) mass-squared — starting from its initial classical (Bunch-Davies) 
value. This means that linear perturbations do not grow in this model. Thus, the model is 
stable, (iii) The results of quantum field theory and Starobinsky's stochastic technique are 
in perfect agreement. 

The effect can be understood physically as follows. Heisenberg's energy-time uncertainty 
principle implies that virtual particle-antiparticle pairs continually emerge from the vacuum 
and then disappear back into it. However, massless particles which are also conformally 
noninvariant have a certain amplitude for appearing with wavelength greater than the inverse 
of the Hubble parameter. In an inflating universe, these virtual pairs are pulled apart by the 
Hubble flow before they find time to annihilate each other. Hence, they become real, recalling 
the analogy with the Hawking radiation. Therefore, one gains particles out of nothing which 
also means that the scalar field strength grows. In fact (Q\^ 2 \Q) = gift + 0(A) m our 
model. Thus, inflationary particle production causes the scalar to undergo a random walk 
such that its average distance from the minimum of the quartic potential increases. This 
makes the vacuum energy larger; hence p > 0. Now recall the covariant stress-energy 
conservation law p = —3H(p + p). Because inflationary particle production causes p > 0, we 
must have p + p < to satisfy this law. Hence, the Weak Energy Condition is violated and 
the equation of state parameter w < — 1. Will this effect be terminated? If the growing of 
the field, which generates the effect, stops, then the effect terminates. There are two causes, 
in the interacting theory, which would yield the growing of the field come to a halt eventually. 
The first cause is the classical restoring force — |y9 3 . This force pushes the field back down 
to the configuration where the potential is minimum, i.e., to ip — 0, as the scalar rises up 
its potential. The second cause is the curvature associated with field being away from the 
minimum of the potential, i.e., |y9 2 which acts like a "mass-squared." Because |(f2|y9 2 |fi) 
is a growing positive real number, the scalar develops a growing positive mass. That cuts 
off particle production, since inflationary particle production requires effective masslessness 
(with respect to the Hubble parameter). Because of these two causes, the field can not 
continue to roll up its position. It must eventually come to a halt. In fact, Starobinsky and 
Yokoyama showed [7] that (Q\ip 2 \fl) does asymptote to a positive constant. 
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APPENDIX A: INTEGRATING THE TWO-LOOP TERMS 



In this appendix we calculate the integral 



- / drf d 3 x'Ml(x;J)u(ri',k)e- S( - s - s! ' ) 



(Al) 

which is a part of Eq. (40). M.%(x;x') is defined in Eq. (32). In the late time regime of 
physical interest, we replace u(rj' : k) (Eq. (15)) with its constant limit u(0,k). Moreover, as 
can be seen below, when ++ and H — terms are added, factors of the Heaviside function 
©(At? — Ar) arise. They require \\x — x'W < 77 — 77', ensuring casuality in the Schwinger- 
Keldish formalism. In the late time limit this means that the spatial plane wave factor 
e ifc-(x-x') can a j go kg d r0 pp ec i j n Eq_ (Al), since 77 — 77' — > in that regime. We, therefore, 
break the integral (Al) into a sum of six terms of the general form Yln=i In(v)i with 

4(77) = -u(0, k)J drf J d 3 x'Ml n (x; x') . (A2) 
Here, we define the integrands 

-ln(/i 2 A^ + ) ln(/i 2 A^_)- 



Ml^x-x' 
M\ 2 {x\x' 
M\ z {x-x' 

Ml A (x;x' 

M 2 2:5 (x-x' 

M 2 2fi (x;x' 
H 2 

2 7 3 2 7T 4 

H 2 A 



2 12 3tt 6 
1H 2 

1H 2 



aa'd 4 
aa') 2 d 2 

'\2a2 



Ax 2 ++ 



Ax 2 



In 



Hel\ fin (fx 2 Ax 2 ++ ) In (fx 2 Ax 2 + _) 



VYu) 



2 n 7T 6 



aa'fd 2 



2fi J V Aa;2 + A4 
In 2 (/i 2 Aa; 2 + ) In 2 (/i 2 A:r 2 _) ' 



A* 2 + 



Ax 2 



2% 6 



/\3 



In 2 (^-H 2 Ax 2 ++ ^ In 2 (^# 2 A:r 2 _) 



Ax 2 



Ax 2 



2 10 3tt 6 



2 9 3 

.3 



ln 3 (^// 2 Ax 2 + ) -\n*(^H 2 Ax 2 + _) 
l — a 2 jln(a)<9 2 - (2 \n(a) + l^Had ^5\x-x') 

39 + 27 ln^)-2vr 2 ] ln(a) ^5 4 (x-x') 



(A3) 
(A4) 
(A5) 

(A6) 
(A7) 



a 4 <! 41n 3 (a) + — ln 2 (a)- 
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y 



+ 



n + 5 



-a 



-(n+l) 



-(n+2) 



- 4 E 



a 



-(n+3) 



2 7 tt 4 ^ 1 81 ^(77 + l) 3 ~ ' ^ (n + 2) 3 ' ^ n(n + 3) 3 ' ^ ) ' (A8) 

v n=l v ' n=l v ' rz=l 

We evaluate ^1(77) and ^3(77) explicitly to illustrate the relevant calculation techniques and 
give only the final answers for the remaining four that can be obtained similarly. The first 
integral is 

n r 

(A9) 



h(rj) = -u(0,k) f drf j d 3 x'Ml 1 (x;x') . 
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It is useful to break up the logarithms in M. 2 .^ as 

hVA* 2 ) = ln(^-) +2 In (|) . (A10) 
We then partially integrate the inverse powers of Ax 2 , using the identities 

:f \n(H-\r-) . (A12) 



Ax 2 4 

Because the integral is over x /At , the derivatives with respect to x^ can be taken outside 
the integral. The remaining integrand possesses only logarithmic singularities. There is no 
distinction between the ++ and H — terms at this stage. We define the temporal and spatial 
intervals At] = rj — rf and Ar = \\x — x'\\, respectively. The ++ and H — terms cancel for 
Arj < 0, so we can restrict the integration to Arj > 0. Then, the logarithms can be expanded 
as 

ln(^-A4 ± ) = ln(^(A?7 2 -Ar 2 )) ± t^Q(Arf- Ar 2 ) . (A13) 
Making the change of variables r = x — x' and performing the angular integrals yield 
IM =-^^«(0,A;)a9 6 £^V^ A Vr 2 [ln(^(Ar ? 2 -Ar 2 ))+21n(^)-l" , 

where the initial time rji = —if -1 . Next, we make the change of variables r = Arj z and 
perform the integration over z, using 

-\- < A14 > 

The result is 

Hv) =~l^^-A^ k ) ad l /V^ 3 [Ar ? 3 (ln(iJAr ? )+ln(^)-^)] . (A15) 

Note that, owing to the factor Ar/ 3 , three of the external derivatives were brought inside 
the integral. This cubic derivative gives 6\n(2/iAr]), when it acts upon the terms inside the 
square bracket. At this stage, one makes the change of variables rj' = —(Ha')' 1 and looks 
up the relevant integral, Eq. (Bl), from Appendix B. Acting on the remaining derivatives 
using do = Ha 2 J^, one obtains 

h(v) = 2^«(0, k)a4 ln(a) - m(f ) + \ + £ ( " ~ ffi ~ 2) aA . (A16) 

^ n=l ~" ' 
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Evaluation of hiv), denned by Eqs. (A2) and (A4), is similar to that of hirf). Using 
Eqs. (A10)-(A14) and (B2) yields 

H 2 , (He\\ , s , f , _ , /2/i\ ^ . /. 1' 



T M = K^H ' fc ) a4 { ln W - ln (£) + 1 ~ ln ( 1_ a) ~ ° J " (A17) 
Next, we explicitly evaluate 

13(77) = -«(0, k) f drj' j d 3 x'Ml :3 (x;x') , (A18) 
where A4 2 3 * s gi yen m Eq. (A5). We break up the logarithm squared in M.\ 3 as 

W = +41,(1) h (*^) + 41n»(|) . (A19) 

Then, we use the identity 

ln 2 (iJ 2 A:r 2 ) d 2 , , , rr2 A 2 , <9 2 , 2 , rr2 A 2 , 2 , A . 

As* = 12 M ^ ^ - T ln + A^ (A20) 

and Eqs. (A12), (A13). To evaluate the radial integral, we make the change of variables 

r = Arjz and use Eq. (A14) and the integral 

■l-z 2 \ 104 7T 2 



J dzz 2 \n 2 (- 



We find 

^3(77) 



H 2 



2 8 3tt 4 



27 9 



M(0,fc)a 2 a / d^a'^g A?f ln'(#Ar7) + Aln(#A77) + 5 



(A21) 



(A22) 



where A = 2 ln (2/////) -11/3 and B = In 2 (2/V#) - (11/3) ln (2fi/H) + 85/18 - tt 2 /6. The 
cubic derivative in the integrand yields 6 ln 2 (2yuA^) — 7r 2 , when it acts upon the terms inside 
the square bracket. Making the change of variables rj' = —(Ha'Y 1 and using Eqs. (B2)-(B6) 
and do = Ha 2 -§-, one obtains 



!M = ^(°> k)a4- l -\n\a)+ [ln(^ ) -l] ln(a) - [± ln(| ) -l + a^] ln( 



ln 



2(1 

~H 



+a 



Infl-^-^l^ri-iV^^, 



a 



a 



n=l 



2/_f\ 7T 2 

IF/ 12 



(A23) 



Evaluation of 14(77), defined by Eqs. (A2) and (A6), is similar to that of h(rj). Using 
Eqs. (A12)-(A14), (A19)-(A21), (B2), (B3), (B6) and (B7) yields 



H 2 



+ 



2 6 7T 4 

l-2a~ 1 + a" 2 



1 7T 2 



h(v) = 7^J«(0, k)a* <{ In 2 (a) - - ln(a) + - + - + 



2 + °" 



In 2 1 





rl3 


7T 2 - 


a 1 + 


-25 




I 8 








3 . 


.16 



16 6 J 



n z 



n=l 



(A24) 
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To evaluate 75(77), we expand In 3 ^H 2 Ax 2 /ij = \n^H 2 Ax 2 /4^ + (1/2) 
Eqs. (A13), (A14), (A21) and (B5)-(B8). The result is 



and use 



+ 



^^u(0,k)a 4 L 3 (a)- 9 -\n 2 (a) + 



15 7T 2 1 , . . 2035 7T 2 

ln(a)-^^ + 



16 



288 6 



71 25 35 _ 2 13 „ 

a H a a 

12 2 4 6 



lnfl 

a 



— -9a _1 + -a _2 -a" 3 
2 2 



ln 2 (l-- 

a 



+ 



671 3tt 2 i 
"48" 2 



! [883 3tt 2 



96 



a- 2 + 



329 7T 2 



144 6 



n=l 



where the Digamma function ?/>(n) = —7 + X)fc=i ^ 1 an< ^ Euler's gamma 7 ~ 0.577. It is 
straightforward to show that the remaining integral 



^(.) = ^(0^)a 4 {ln 3 (a) + f ln 2 (a) + [^ln(|)-f + ^] ln(a)-l ( 



» V- n + 5 (n+1) U L __ U__ 

^^(n + 1) 3 ^(n + 2) 3 ^n(n + 3) 

n=l v ' n=l v ' n=l 

Summing the six terms gives the total two-loop contribution in Eq. (40) 



-0+2) 
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27, /2/i\ 189 7T 2 
'77) T 



In (a) 



+ 



8 V if/ 8 \HJ 



'2fi\ 205 13 2 

+ T^ 7r + 



13 a" 



144 48 
79 _., 13 



115 25 



a _1 +a -2 -— a -3 

48 4 12 



In 



, 2/-, A ' y -1 lo _o 

ln V-a) + Ts a ~T2 a + 
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n=l v ' 



Q ~ fl -(n+2) 9 x 
"9 2^ („M0\3 9 



a -(n+3) 3 ~ 
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APPENDIX B: USEFUL INTEGRAL IDENTITIES 



In Appendix A, to calculate the temporal integrations over 77', we make the change of 
variables rf = —(if a') -1 , and use the following integral identities: 

(1-0 



A a' \a' aJ 2 w ^ 

n=l 

^ln(i-i) =-aMa) + (a-l)ln(l-i), 



(Bl) 
(B2) 



24 



r daVin G4) = 4 in w4 ( ° 2 - a)+ 5( a2 - i ) in ( i 4) • 

f , , „, /l 1\ a 3 , , . a 3 a 2 a 1. „ <N , A. 1 



! a' Va' a 



2 

ln 3 (a) + ^ln(a) + 2^ 



n=l 



n z 



■0(n)+7- 



nJ 



(B3) 
(B4) 
(B5) 



where ip{n) is the Digamma function and 7 is the Euler's constant, as defined in Appendix A, 



(1 --P + (-Dln 2 (l-i 



(B6) 



rfa'ln 2 ^-^) = aln 2 (a) + 2a^ 
f 1 da'a'ln 2 (^--) = ^ ln 2 (a) + a 2 ln(a) +a 2 ^ " } - a(a-l) ln(l--) 

n=l 

4 (a2 - i)ln2 ( i 4)- (B7 » 

1/1 n=l 

-^(a-l)(3a+l)ln(l-i)+i(a 3 -l)ln 2 (l-i) . (B8) 
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